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We investigate the Bianchi I cosmology with the homogeneous SU(2) Yang-Mills field governed by 
the non-Abelian Born-Infeld action. A similar system with the standard Einstein- Yang-Mills (EYM) 
action is known to exhibit chaotic behavior induced by the Yang-Mills field. When the action is 
replaced by the Born-Infeld- type non-Abehan action (NBI), the chaos-order transition is observed 
in the high energy region. This is interpreted as a smothering effect due to (non-perturbative in a') 
string corrections to the classical EYM action. We give a numerical evidence for the chaos-order 
transition and present an analytical proof of regularity of color oscillations in the limit of strong 
Born-Infeld non-hnearity. We also perform some general analysis of the Bianchi I NBI cosmology 
and derive an exact solution in the case when only the U(l) component of the Yang-Mills field is 
excited. Our new exact solution generalizes the Rosen solution of the Bianchi I Einstein-Maxwell 
cosmology to the U(l) Einstein-Born-Infeld theory. 
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I. INTRODUCTION 

One of the key questions in theoretical cosmology is 
whether the space-time metric near the singularity is 
regular or chaotic. As was shown by Bclinskii, Kha- 
latnikov and Lifshitz (BKL), the generic solution of the 
four-dimensional vacuum Einstein equations exhibits an 
oscillating behavior Jl| which was later qualified as essen- 
tially chaotic (see and references therein). Recently 
the issue of chaos in the early universe received a re- 
newed attention due to discovery that the antisymmetric 
form fields in ten and eleven-dimensional supergravities 
imply chaos 0, 0, 1^ lU ■ Namely, it was shown that 
the general solution near a space-like singularity of the 
Einstein-dilaton-p-form field equations exhibits an oscil- 
latory behavior of the BKL type. However the issue of 
chaos in superstring cosmology is not completely solved 
yet, since these considerations are based on the lowest 
in a' level of the string theory. To go beyond this ap- 
proximation in the closed string theories is difficult, since 
no exact in a' effective action is known. Moreover, the 
quadratic curvature corrections obtained perturbatively 
are likely to disfavor chaos M, 1^ . Recently damping of 
the BKL oscillations was observed also within the brane- 
world scenario p^ . 

The role of non-perturbative string corrections is eas- 
ier to study within the open string context (or D-branes) 
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where an exact in a' effective action is known in the slow 
varying field approximation under the form of the Born- 
Infeld (BI) action [Hill- A non-Abelian generalization 
of the Born-Infeld action was suggested in ref. 13] . Con- 
sequently, one can explore whether the Yang-Mills (YM) 
chaos extensively studied for homogeneous fields in fiat 
space-time and in Bianchi I cosmology persists when the 
usual YM action is replaced by the non-Abelian Born- 
Infeld action (NBI). This would probe the effect of the 
string non-locality on the issue of chaos. 

Classical YM fields governed by the ordinary quadratic 
action exhibit chaotic behavior in various situations 
[11 m. The simplest case is that of the homogeneous 
YM fields depending only on time in the fiat space-time 
[TgL IT^ IT^ : when only two YM components are ex- 
cited, the problem is reduced to the well-known two- 
dimensional hyperbolic system H — (p^ + + 2;^y^)/2 
which is chaotic. Furthermore, in the lattice simulations 
of the inhomogeneous YM system, one observes the en- 
ergy flow from the infrared to the ultraviolet region |l9j |. 
Therefore, it is believed that the chaotic behavior is typ- 
ical for the purely classical YM equations, one of the 
arguments being the absence of solitons in this theory. 
In addition, it is known that adding the Higgs field to 
the YM theory leads to stabilization of chaos in the ho- 
mogeneous systems. In this case the hyperbolic model is 
replaced by the system of coupled harmonic oscillators 
which is regular in the weak coupling regime. 

In the context of YM fields in the presence of grav- 
ity, it is of interest to pint the following. The YM 
field has violent oscillating behavior near the singular- 
ity of the Einstein- Yang-Mills (EYM) black holes '2^, 
but the oscillations are not chaotic. In the domain of 
cosmology, some homogeneous models, such as an ax- 
isymmetric Bianchi I, the YM chaos unambiguously per- 
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sists j2l|, |22, 123 , though in principle gravity leads to a 
smoothing of the chaotic behavior. On the other hand, 
the Born-Infeld effect on the flat-space dynamics of the 
homogeneous axisymmetric YM field was shown to pro- 
vide a chaos-order transition [2^, so it can be expected 
that in the gravity coupled case this effect will be even 
more pronounced. Note that in the related investigation 
of the behavior of the YM field inside black holes [2^, 
it was found that violent YM oscillations disappear once 
the quadratic YM action is replaced by the NBI action. 

In this paper, we study in detail the axisymmetric 
Bianchi I cosmology with the YM field governed by 
the NBI action. Though the system becomes much 
more complicated when the ordinary YM Lagrangian 
is replaced by the non-Abelian Born-Infeld (NBI) La- 
grangian, we show that the equations can be considerably 
simplified in some physically interesting limiting cases, 
and even admit some exact solutions. The generic solu- 
tion exhibits a transition from the YM chaos to the reg- 
ular oscillating regime when moving backward in time. 



II. GENERAL SETTING 



As was discussed recently [l3;l2fll23,l2S|, the definition 
of the NBI action is ambiguous. One can start with the 
U(I) BI action presented either in the determinant form 



or in the equivalent (in four space-time dimensions) 
"square root" form 



S = 



1 

l"67r 




2/32 



16^4 



^gd^x. (2) 



In the non-Abelian case of the matrix- valued i^^jy, the 
trace over gauge matrices must be specified. One par- 
ticular definition is due to Tseytlin ^T"]. A symmetrized 
trace was therein introduced, prescribing a symmetriza- 
tion of all products of F^i, in the power expansion of 
the determinant ^ before the trace is taken. Inside 
the symmetrized series expansion the gauge generators 
effectively commute, so both the determinant and 
the square root ((SJ forms are equivalent. This property 
does not hold for other trace prescriptions, e.g., an ordi- 
nary trace. In the latter case it is common to apply the 
trace to the square root form 0. Note that string the- 
ory seems to require the symmetrized trace definition in 
the lower orders of the perturbation theory |23, , 
while higher order corrections seem to violate this pre- 
scription [33, 1^ 113, ■ Here we choose the "square 
root/ordinary trace" Lagrangian just for its simplicity. It 
is worth noting that in the static case discussed rece ntly 
both in the ordinary 0| and the symmetrized trace [23 
versions, qualitative features of solutions turn out to be 



the same. Thus we choose the action of the Einstein-NBI 
system in the following form 



S 



I 

47r 



1 

AG 



-gd'^x, (3) 



where R is the scalar curvature, /3 is the BI critical field 
strength and 



n 



1 



pa pMi' , 
2/32 " 



1 



16/3' 



iF^.Fry 



(4) 



The limit /3 — > 00 corresponds to the standard EYM 
theory with the action 



-gd'^x. (5) 



We consider an axially symmetric Bianchi I space-time 
described by the line element 



N^dt^ - b^{dx^ + dy^) - c^dz"^ 



(6) 



where functions TV, b and c depend on time t. In the 
YM case this problem was studied previously by Parian 
and Kunzle [21I |2^ and Barrow and Levin j23|. The 

gauge field compatible with the space-time symmetry is 
parameterized by two functions u, v of time 



A = Tiudx + T2udy -\- T^vdz, 



(7) 



where SU(2) generators are normalized according to 
[ri,T2] = iT^. The corresponding field strength matrix- 
valued two-form is 

F = u{Tidt hdx + T2dt A dy) + uTgdt A dz 
+ u'^Tsdx Ady + uv{T2dz A dx + Tidy A dz). (8) 

Integrating over the 3-space, we obtain the following one 
dimensional Lagrangian 



L = - 



1 b{hc + 2cb) 

2 GN 



0^Nb'^c{n~\), 



where now 



7^ = 




Nb^c 



(9) 

(10) 

(11) 
(12) 



The quantity J-' is the YM Lagrangian, and it is conve- 
nient to present it as a difference of kinetic and potential 
terms 



J" = T-U, 
2u2 

T = 
U = 



+ 



iV262 
2'u2w2 



62 c2 



u 



(13) 
(14) 

(15) 
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Note that from two coupling parameters entering the ac- 
tion, G and /?, one can be eHminated by an appropriate 
rescahng. In what follows we set G = 1. 

The Einstein equations can be derived by variation of 
the one-dimensional action over N, b, c. Variation over N 
gives the Hamiltonian constraint 



V- -n 



(16) 



where Ti. reads in the synchronous gauge A'^ = 1: 

n^hibc + 2cb) + ^[p\n-i)-u]. (17) 



Fixing this gauge from now on, we obtain the remaining 
Einstein equations: 



b be c 
b be c 
b 162 
6 ^ 262 



2/32(7^-l) 



n \b^ b^c^ 54 /?2 

,-,2 ,,2, ,2 g2 



1 

n 



2 2 

U V 

'6V 



1(19) 



The equations for the YM field can be presented in the 
following form 



n d 
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( . uvG 


c dt 
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' 

cJl— 

dt 


cTZ 


( cuvQ 
r /32 



(iiv + vu)Q V? uv^ 



cfP 52 c2 



2cuuQ 



?o) 

(21) 



The energy-momentum tensor has the following com- 
ponents. The energy density is given by 



/32 



2^2p2 



•^0 



47^7^ 



47r' 



(22) 



the pressure in the plane orthogonal to the symmetry 
axis is 



-r2*2_^2 p2 



and the pressure along the axis of the symmetry is 
2n* - - /32 



Vz 



47^7^ 



47r' 



(24) 



the Lagrangian remains invariant. Moreover, under a 
separate rescaling in the 6, u sector 



A&, u — !■ Aw, 



(25) 



the Lagrangian scales as A2, and under the transforma- 
tion 



Ac, V Aw, 



(26) 



as A. The corresponding reduction of the EYM system 
is achieved by an introduction of new variables invariant 
under the above rescalings. Following Barrow and Levin 
|23| . whose notation we will adopt in what follows (note 
that in Ref. |2^ another convention SttG = 1 is used), 
we introduce the volume and shear variables 



fb\ 



together with the associated Hubble parameters 

X 



X 



as well as the scaled Yang-Mills variables 

* = r = -. 

c 

It is also convenient to use the scaled derivatives 



= 7, Hr = -, 
c 



which are related to ^, F via 

n^ = i' + {Ha + H^) 

Hr = f + (i/a - 2H^) F. 



(27) 



(28) 



(29) 



(30) 



(31) 
(32) 



Note that these are not the momenta conjugate to ^,F, 
the corresponding canonical momenta being 



2a3 / ^-F^A 



(33) 
(34) 



In terms of the new variables the Hamiltonian con- 
straint reads 



III. REDUCTION OF ORDER 

The above system of equations look as a dynami- 
cal system of the eight-order in the presence of a con- 
straint. However, it possesses additional scaling symme- 
tries which can be used to reduce the system order by 
two (for the EYM action this possibility was noticed by 
Darian and Kunzle '2l'|). It is easy to check that under 
a scaling transformation 



Afo, 



A" 



Au, 



A-2„^ 



-(i^2-i^2)+^2_^^2_^^2(-^2_^2F2)]7^"l =0. (35) 

The functions T, U and Q entering TZ now take the form 



T = 2H^+Hf, 
U = ^'* + 2*2r2, 

g = *(2H3,F + Hr*) 



(36) 
(37) 
(38) 



From the Einstein equations one can derive two first 
order equations for the Hubble parameters, which are lin- 
ear in derivatives. Taking the sum of the Eqs. (|18|l . (|19|) 
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and the constraint equation (|17|l , one obtains the follow- 
ing simple equation for Ha'. 



and bring the Einstein equations into the form 



Ha + iHl + —{T + U) =0. 



(39) 



Similarly, taking twice the second Einstein equation H19|l 
and subtracting (|18(l we get 



Hb + HtiHh-H,) = --(n| + *2r2), 



(42) 
(43) 



3i/x^a- 3^(n2 -n| + *4-*2r2) = o. (40) 



Thus the Einstein equations reduce to two first order 
equations in the presence of a constraint. 

Alternatively, one can introduce the Hubble factors 
with respect to b and c: 



H, 



(41) 



with the Hamiltonian constraint 



^HtiH, + - i [p^l -TZ) + U] 



(44) 



In addition, we have two second order equations for 
the YM fields which read in terms of the new variables 



di 



1 



d 
di 



n 



nr4 



1 



2 

n 



g(n.i,r + nr^) 

/32 



0, 
0. 



(45) 
(46) 



IV. YM LIMIT 

In the YM limit (3 ^ oo the square-root factor in the 
above formulas should be replaced according to the rela- 
tion 

1 



lim I3\n 



1 



(47) 



The main qualitative difference between EYM and ENBI 
theories lies in the fact that the standard YM action is 
scale-invariant (though not the EYM one) contrary to the 
NBI case. This leads to a partial decoupling of the YM 
dynamics from that of the space-time. Given eq. (|47|l . 
the constraint equation simplifies to 

^[3{H'a-Hl)-{T + U)]=0. (48) 

Combining this with (|39|l . one finds that one of the Ein- 
stein equations fully decouples and reduces to the vacuum 
form: 



H„ + Hi 



2Hf. 



0. 



(49) 



However, the shear remains coupled to matter and obeys 
the equation 

,4 



H-^ + 3H^Ha + Ha — H^ — H|, 
Finally, the YM field equations become 



tl^ + {Hb+Hc)Il^+^{-^^ + T^) = 0, 
Hr + 2iJbHr + 2^'2r = 0, 



(50) 

(51) 
(52) 



where the definitions 1)31(1 . 1(3 2|l have to be used. 

The Hamiltonian form of the EYM equations can be 
further simplified using an exponential parametrization 
of the volume and shear variables 



a = e , x = e'. 
The canonical momenta conjugate to a, 7 are 

Pa = -3e3"d, = 3e3"7, 
while the YM momenta ((33(1 . ((34(1 simplify to 



(53) 
(54) 
(55) 



P^, = 2a^U^, Pr=a'Ur. 

The Hamiltonian constraint l(17() for the EYM system in 
terms of the momentum variables reads 



n 



-3a 



6 ^ " 



where the potential is given by the Eq. ((37(1 . 



(56) 



V. U(l) CASE 

Consider the special case when only the w-component 
of the YM field is excited, corresponding to the U(l) 
subgroup of the gauge group. The Einstein equations 
(1321, 63 reduce to 



Hh 



H, 



Hb{Hb 

Hc{Hb 



Hr 



Hr 



0, 



n 



rH 



r, 



(57) 
(58) 
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and the Hamiltonian constraint is 

Hb{Hb + 2Ha) ^ 2/32 (7^-l - 1) 
Integrating the BI field equation 



one obtains 



dt \ n 



n 



0, 



(59) 



(60) 



(61) 



where bo is an integration constant, so that 

1 2h 



/32 VTT^ 



' (62) 



It is easy to see that the Einstein equation H57|) is equiv- 
alent to 



b e 

which immediately gives a relation 
b — kc. 



(63) 



(64) 



where fc is a second integration constant. Now the con- 
straint equation becomes the following separated equa- 
tion for the function b{t): 



(65) 



while the second Einstein equation (|58|l is its time deriva- 
tive. The right hand side of this equation is positively 
definite. It follows that the system has no bounces. In- 
deed, if Hb = 0, from the Eq. (|^ it follows that Hb = 0, 
which contradicts the Eq. (|65|l . 

We can solve the Eq. (|65(l considering instead of b{t) 
an inverse function t{b). Then 



(66) 



where t' = dt/db. The equation for t{b) following from 
reads 



This is the linear first order equation for the function 

<b) = (1A')^ (68) 

namely. 



0. (69) 



Its solution reads 

2/32 




i+3.~i\b^db+'-^^ 



/32 64 



(70) 



where bi is a third integration constant. An integration 
can be done in terms of the hypergeometric function [s^ : 



2/32 



/64 



/32 

8f3bbQ 



2u2 



F 
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3' 4' 4' 1 + x2 



bi 



(71) 



where bi ^ 61 is another constant. Now, according to 
the inverse function to the required solution is given 
by the integral 



tib) 



db 



+ to, 



(72) 



where to is the last integration constant in t his p rocess. 
Our solution generalizes the Rosen solution ^8] to the 
Einstein-Born-Infeld theory. 

Near the singularity z « bi/b, so one has 



Integrating Eq. H72|l one obtains 

& = (M)'/^ 

and then from Eq. 164(1 



2b 



2/3 



3k 



-t 



-1/3 



(73) 



(74) 



(75) 



Hence, we obtain a cigar singularity. 

In the Maxwell case the situation is different. Indeed, 
in the limit [3 ^ 00 one has 



4fe| 61 
b' 



(76) 



Since z should remain positive, the region of b is limited 
by 



b > 6min = -r^- 

bi 

Combining the Eqs. JT^J), we obtain 



46mi: 



bit 



(77) 



(78) 



This is a pancake singularity. Thus, the BI non-linearity 
modifies the singularity from a pancake to a cigar type. 



6 



VI. SINGULARITY STRUCTURE 

Consider now the general solution near the cosmo- 
logical singularity. It turns out that except for a spe- 
cial isotropic solution b = c = a, previously studied in 
HE HE QE ^1 1 generic solutions have the same metric 
singularities as the vacuum Bianchi I solutions. Near the 
pancake singularity the solution is not analytic in terms 
of t, but in terms of t^^^. In fact, one finds the follow- 
ing Laurent expansion containing four free parameters 
p,q,r,s: 



1 2rq + sp 2 



3t 



9pr 



2/3 



(79) 



2/3 



0(1), 
0{t). 



(81) 
(82) 



The F-component of the YM field vanishes at i = 0, while 
is singular. The scale factor a and the shear x near 

the pancake singularity both behave as 0(t^/^). 

Near the cigar singularity the solution has a Laurent 

expansion in terms of t: 



1 



= 



37^l 

2f2p3 



\-0{t), (83) 
■ + 0{t), (84) 
t + 0{t^), (85) 



ft'^ +S + 



7^l 



0(i), 



(86) 



where the quantity TZi is the leading term in an expansion 
of the NBI square root: 



/?2 



(87) 

The scale factor and the shear have the following expan- 
sions: 

a = ai {t'/' + ^^1^ t'^' + 0{ty^)^ , (88) 

X = XI (^t-V3 + ?!!|l±l!i2/3 + 0(i^/3)^ _ (89) 

The quantities p, q, r, s (p, g, f, s) are independent 
free parameters which, together with an arbitrariness as- 
sociated with a time shift, provide five constants needed 
to specify the generic solution for both singularity types. 

VII. SOLUTION IN THE LIMIT /3 = 

In order to better understand the effect of the BI non- 
linearities on the gauge field dynamics let us first study 



the strong field limit F ^ /3, or, formally, f3 0. The 
leading term in the square root (|10|l containing the pseu- 
doscalar invariant Q is negative definite. Therefore, im- 
posing the square root TZ to be real-valued in the limit 
(3 —+ may be ensured only if Q tends to zero, in which 
case 



*nr + 2rn* = o. 



(90) 



One can show that this condition is compatible indeed 
with the equations of motion as /3 — > 0. 

Given the condition (|9()|l . the square root term will 
read 



V^^ + 2r^vi/^-n^-2n^ ^ ^^^^ 

The right hand sides of the Einstein equations (|39|l . H40() 
tend to zero, so the gravitational degrees of freedom de- 
couple 



(92) 



and the gravitational constraint assumes the vacuum 
form as well 



Hi 



0. 



(93) 



Decoupling of gravity means that in the limit /3 — *■ 
the metric is given by the vacuum Kasner solution either 
of a cigar type 



Ha — — 



1 

3t' 



or a pancake type 



^ it 



(94) 



(95) 



where we set the singularity dXt — Q. 

Substituting the explicit expressions for the Hubble 
and shear parameters one finds another constraint 



C = const. 



(96) 



and thus in the remaining equations one can express all 
the gauge field variables either in terms of F, Ilr, or in 
terms of 11^. One simple consequence of this con- 
straint is that in the non-trivial case C ^ the variables 
and F can not have zeroes except for the singularity, 
and thus should preserve their signs. From the NBI field 
equations one then finds 



cigar. 



2* 



= — — , pancake. 



(97) 
(98) 



In both cases the dynamical equation for 'J will be of the 
form 



4' = /(*,*,C,t) 



(99) 
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with some function /. It describes oscillations with a de- 
creasing amplitude. The second YM variable T is related 
to 5* algebraically via constraints H90|l . H96|l and there- 
fore oscillates with the same frequency exactly in an an- 
tiphase. Oscillations are fully regular, so no YM chaos 
can persist in the regime of the strong BI non-linearity. 

The general solution near the pancake singularity can 
be expanded with respect to the variable r = t^^^: 

3C 



r = plT + Vecpi + — + qiT^ + oix^), (100) 



* = 



c 



^ r-i + ^^ + O(rUl01) 



2pi 



where pi and qi are free parameters. 

Near the cigar singularity the solution can be expanded 
in terms of t: 

r = PI t-i + q, + ^A_Pl£. t + o{t^), (102) 

6pi 



C 



Cq 



^-t-rOif). (103) 
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VIII. CHAOS-ORDER TRANSITION 

Now we address the problem numerically. Various 
methods were suggested to study a chaotic behavior in 
the context of gravity, where the absence of the canoni- 
cal time variable prevents a straightforward use of such 
convenient tools as the Lyapunov exponents (however, 
see 01)- cds,% of the conformally invariant YM 

Lagrangian, one can use the approach of Ref. [2^ to 
separate the dynamics of the YM field from the gravita- 
tional expansion and then apply the invariant technique 
of chaotic scattering. For systems exhibiting chaotic be- 
havior the set of all periodic orbits has fractal structure 
invariant under coordinate reparameterizations. 

However these methods become problematic in our 
case, where the conformal invariance is absent from 
the matter action. Apart from the special asymptotic 
regimes, it is not possible to separate the YM dynamics 
from the metric evolution. Though it can be done for 
the high YM intensity, it turns out that the time interval 
in the actual evolution where this regime holds is rather 
small for /3 of the order of unity or greater. The YM vari- 
ables perform only a small number of oscillations during 
the epoch of high field intensity, the field then being fast 
diluted. However, if we set the parameter /? sufficiently 
small, the time spent by the system in the highly non- 
linear region will be large enough, and in this case the 
chaos-order transition is unambiguously manifest. 

A numerical analysis of the system for small values of 
/? reveals the following. While the gauge field strength is 
considerably greater then the critical field both condi- 
tions (|90|) and H9()|) approximately hold, and the dynamics 
of the gauge field qualitatively coincides with that dis- 
cussed in the previous section for /3 = 0. Both variables 



and r perform nonlinear oscillations with decreasing 
amplitude in an antiphase with respect to each other and 
without crossing zero. In this region the dynamics is fully 
regular, while the evolution of the metric is mainly gov- 
erned by purely vacuum terms and is close to the vacuum 
Kasner solution. 

To test the validity of an approximate description of 
the system in terms of the limiting /? = solution we 
check the constraint equations H90I96I) in the case of small 
but finite /?. The first constraint is the condition of 
smallness of the pseudoscalar YM invariant ^ (?T . 
The second constraint was directly checked numerically. 
Fig- n illustrates the situation for the cigar solution with 
/3 — 10^^. From this figure one can see that the right 
hand side of Eq. H9()|l evolves on time scales much larger 
than the period of oscillations of the gauge field vari- 
ables. Under the overall volume expansion, the YM en- 
ergy density falls down and the role of the BI nonlin- 
earity decreases. At the same time the matter terms in 
the Einstein equation become more significant forcing, 
in particular, to decrease the shear anisotropy much 
faster than the Hubble parameter Ha- 

These features are illustrated in Figs.|2[El Fig.|21shows 

the early regular evolution for /? = 2 • 10""^ and the cigar- 
type singularity. Both variables ^' and F oscillate in the 
positive region. The behavior of the shear anisotropy 77^ 
is smooth. Fig. O demonstrates the same solution at late 
time. One can see that the dynamics of the gauge field 
becomes essentially chaotic. The function coupled to 
matter performs chaotic oscillations with decreasing (as 
compared to the Hubble parameter Hg) amplitude. The 
first zeroes of the gauge functions ^E* and F serve as an 
approximate boundary separating highly nonlinear evo- 
lution from the region of the chaotic regime. The actual 
type of singularity (pancake or cigar), plays a relatively 
small role in both the regular and strongly chaotic phases 
except for the small vicinity of the singularity. This is il- 
lustrated in Fig.2|which shows the solution with pancake 
singularity obtained from the solution shown in Fig.|2by 
changing the sign of with initial conditions which 
were set at t = 10. 

The Hubble parameter Ha does not exhibit chaotic be- 
havior. It can be presented as Ha = h{t)t~^, with some 
slowly varying smooth function h{t). Numerical curves 
h(t) are shown in Fig. |31 for various (3. This function 
interpolates between the value 1/3 at t = (vacuum 
Kasner solution) and the value 1/2 at t = oo correspond- 
ing to the isotropic "hot universe" cosmology. However, 
for small values of /3, when the system stays in a highly 
nonlinear regime for a considerable time interval, there 
is a region where h{t) is greater than 1/2. This feature 
can be explained using the results of the FRW-BI model 
and was firstly presented in ref. 35| (Subsequent anal- 
ysis can be found in ref. [11101 41].). Eq. implies 
that once the contribution of the anisotropy term H^ de- 
creases (i.e. the solution undergoes isotropisation), the 
Hamiltonian constraint tends to the Friedmann equation. 
In the FRW case [s^ one can derive an equation of state 
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for the NBI matter which interpolates between that for 
conformal matter e = l/3p and the "string fluid" equa- 
tion e = — l/3p in the highly nonlinear regime. The 
latter corresponds to the value h = 1, which is, however, 
never achieved in the anisotropic case which we investi- 
gate here. 

IX. DISCUSSION 

The main goal of this paper was to test the non- 
perturbative effects of superstring theory on the issue 
of chaos in cosmology. At least three different patterns 
of chaotic behavior in cosmology were identified. The 
first is the billiard-type behavior which is manifest in 
the Bianchi IX pure gravity and its supergravity (includ- 
ing multidimensional cases) generalizations. The second 
is the bouncing behavior of the FRW-scalar field cos- 
mology. The third type is the matter-dominated chaos 
of the Bianchi I EYM cosmology, and it is this type of 
chaos which was investigated here (Recently an interest- 
ing analysis was performed 43] of the YM field behavior 
in more general type A Bianchi space-times showing that 
basic features of the YM chaos persist there as well). 
From these three patterns the last one is the most ap- 
propriate for testing the superstring non-locality effects 
accumulated in the BI non-Abelian action. Our results 
clearly demonstrate disappearance of chaos in the high 
energy regime. 

From a mathematical viewpoint, it is worth noting that 
the Einstein-NBI system of equations admits a reduction 
of order due to the presence of scaling symmetries simi- 
larly to the EYM case. Moreover, in the strong BI regime 
the axisymmetric Bianchi I NBI system can be reduced 
further due to existence of two additional asymptotic in- 
tegrals of motion. This limit is characterized by the dy- 
namical vanishing of the pseudoscalar quadratic invariant 
of the YM field. This simplifies dynamics considerably 



and leads to a decoupling of the gravitational degrees of 
freedom. Color oscillations are still governed by the BI 
non-linearity and are reducible to the one- variable second 
order system predicting perfectly regular behavior. 

Numerical experiments shows that the system behavior 
for sufficiently small (3 consists of a regular phase in the 
high energy region near the singularity and the chaotic 
phase at later time. The regular phase is qualitatively 
similar to that described by the P — approximate de- 
scription. The chaos-order transition is observed when 
one is moving backward in time towards the singularity. 
The singularity itself is either of a cigar or a pancake 
type, as in the vacuum Bianchi I case, though the YM 
field does not tend to the vacuum configuration. Thus 
the non-perturbative in a' string corrections to the YM 
action suppress the YM chaos which takes place at lower 
energies where dynamics of the YM field is governed by 
the ordinary quadratic action. 

In the case of only an Abelian component excited, 
an exact analytic solution of the Einstein-BI system 
was found which generalizes the Rosen solution to the 
Einstein-Maxwell equations. It also exhibits a differ- 
ent behavior in the singularity as compared with the 
Einstein-Maxwell case. 
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FIG. 2: The solution for the /3 = 2 • 10 ^, regular phase, a cigar singularity. The solid line — Tt^^^, the dashed Une — 'it^^^, 
the dotted line — H^/Ha- 
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FIG. 4: The solution for /3 = 2 ■ 10 a regular phase, a pancake singularity. The solid line — Ft^^^, the dashed line — 'ift^^^, 
the dotted line H^/ Ha- 




FIG. 5: The behaviour of h{t) = tHa for various values of f3. 



